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I. INTRODUCTION 



In the preceding paper |I| we extended our earlier work describing interactions be- 
tween a two-level system and a shape-invariant system. Here we present another 
generalization. 

The model studied in Ref. is a generalization of the Jaynes-Cummings model ||. 
In the standard Jaynes-Cummings model the "field" is described by a harmonic oscillator. 
In our generalization it can be described by any shape-invariant system. In developing this 
model we made extensive use of the algebraic approach j|,[5| to the supersymmetric quantum 
mechanics J7J. In this paper we further generalize the model to one with intensity-dependent 
interactions. 

The standard Jaynes-Cummings model is an idealized model describing the interaction of 
matter with electromagnetic radiation. A variant of the Jaynes-Cummings model takes the 
coupling between matter and the radiation to depend on the intensity of the electromagnetic 
field P-|IT|. This model has great relevance since this kind of interaction means that the 
coupling is proportional to the amplitude of the field which is a very simple case of a 
nonlinear interaction corresponding to a more realistic physical situation. The results of 
this model can also give insight into the behavior of other quantum systems with strong 
nonlinear interactions. 



II. THE GENERALIZED INTENSITY-DEPENDENT AND NON-RESONANT 

JAYNES-CUMMINGS HAMILTONIAN 

The expression of the intensity-dependent and non- resonant Jaynes-Cummings Hamilto- 
nian can be written as 

H = A^A + ^ [A, i f ] (<t 3 + 1) + a (a+A \I~A\A + a_/S^ + hAa 3 , (2.1) 

where a is a constant related with the coupling strength, A is a constant related with the 
detuning of the system and &i, with i = 1, 2, and 3, are the Pauli matrices. 

However, the harmonic oscillator systems, used in this context, is only the simplest 
example of supersymmetric and shape-invariant potential. Our goal at this point is to 
generalize that Hamiltonian for all supersymmetric and shape-invariant systems. With this 
purpose we introduce the operators 

S = a + A + <7_i f (2.2a) 



S^Mpti+LVitiit, (2.2b) 

where 

°±=\ (^i ± ^2) ' ( 2 - 3 ) 

and, now, the operators A and A^ satisfy the shape invariance condition 0. Using this 
definition we can decompose the Jaynes-Cummings Hamiltonian in the form 
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H = H + H 



int 5 



(2.4) 



where 



H Q = S , 

H int = a Si + hA a 3 . 



(2.5a) 
(2.5b) 



We search for the eigenstates of H and, in this case, it is more convenient to work with its 
5-operator expressions, which can be written as || 



S 2 : 

Hj n i = a 



TB_B + T^ 



" 




H 2 


" 


B+B__ 




. 


Hi. 



p 



B + B_ B+Tt 



B + B. 

P 



a 



P TB_ V#i 

H l B+f t -f3 



(2.6a) 
(2.6b) 



where (3 = HA/ a. We use the same notation as the preceding paper flj. There we show 
that the states 



m I 



T 
±1 



a 



Cg> I m) 
m + 1) 



(±) 

m+l 



m = 0,l,2, ••• 



are the eigenstates of the operator S 2 

I *m / °m+l I / • 



(2.7) 



(2.8) 



In this case C^Jn+i = C^]m+i [R{ a i)i R(o>2), R( a 3), ■ ■ ■] are auxiliary coefficients and, | m) 
and | m + 1) are the abbreviated notation for the states | ?/> m ) and | ip m +i) [|IJ. 

At this point, we observe that the wave-state ortho normalization conditions imply in the 
following relations among the C's real coefficients 



(±) 

m+l 



(#(±) I *(±)\ = [c (±) l + [c 

\ m I m I ^m ' < 

/u/W 1 vi>( ± )\ - rWrW _ r {±) r iT) - n 

\*m I W m / — U m '-'m °m+l U m+l — U 



(2.9a) 
(2.9b) 



Now, considering that S 2 and Hj n4 commute then it is possible to find a common set of 
eigenstates. We can use this fact to determine the eigenvalues of Hj n4 and the relations 
among the C's coefficients. For that we need to calculate 



H 



int 



#(±)\ = A (±) I 

mi m I m / 



where A^ are the eigenvalues to be determined. Using the Eqs. ( |2.2j ), 
last eigenvalue equation can be rewritten in a matrix form as 



and 



a 



P TB_ \JH X 

Hx fi+Tt -p 



T 
±1 



Cg> I m) 



A (±) 



a 



Cg> I m) 
m+l) 



(±) 

m+l 



(2.10) 
the 

(2.11) 



Since the C's coefficients commute with the A or A* operators, then the last matrix equation 
permits to obtain the following equations 
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'<*P - (TC^) f\m)±a C^TB- \[W X | m + 1) 



a (fC^T*) yJH x B + \ m) t [a(3 + A^] \ m + 1) = C 



Introducing the operator [12 



Qt= (B + BY 1/2 B, 



,-1/2 



one can write the normalized eigenstate of Hi as 

| m) = (Qt) m | o) 
and, with Eqs. (|2.13|) and ( |2.14| ) we can show that [ 



B + | m) = J£ m+ i | m + 1) , 



T5_ | m + 1) = JS m+1 T\m). 



Therefore, we have that 



TB_ V^i | rn + 1) = T5_ V^+i | m+1) 



£ m+ i T5_ | m + 1) 



and 



Hi B + | m) = V ifi J£ m+ i | m + 1) 



£ TO+1 V-^i I m + l) 

£ m +l | m + 1) , 



Using Eqs. ( gTg ) and ( gT?p , then Eqs. ( g7Lj| ) take the form 

{ [a(3 - Ag=>] (f C^f f ) ± a £ m+1 C^ijf | m) = 
{a S m+ i (f C&T t) =F [a/3 + A^J C^} | m + 1) = 

From Eqs. (|218|) it follows that 



and 



a 



c±) 



^ +1 + /3 2 T /3' 



m+l 



m+l 



Eqs. Q and (HO) im p!y that 



a 



(±) 

m+l 



C (=F) 



(2.21^ 



and the eigenstates and eigenvalues of the generalized intensity-dependent and non-resonant 
Jaynes-Cummings Hamiltonian can be written as 



E<& = £ m+1 ± J a 2 El 



■"m+l 



h 2 A 2 



and 



c£0 1 m ) 

CW I m+l) 



a) The Intensity-Dependent Resonant Limit 



I ¥±h = 

1 m 1 


T 


" 




.0 


±1_ 





m 



0,1,2, 



(2.22) 



(2.23) 



From these general results we can verify two simple limiting cases. The first one corre- 
sponds to the resonant situation, which is for A = ((3 — 0). Using these conditions into 
Eqs. (|2.20|) and (|2.22 ) and Eqs. ( j2.9| ) we can promptly conclude that 



(2.24) 



and 



r (±) 



f n{±)f\ = M±) = J_ 



(2.25) 



Therefore the intensity-dependent Jaynes-Cummings resonant eigenstate is given by 



m ' V2 



T 
±1 



| m) 
m+l) 



m 



0,1,2, 



(2.26) 



If we compare this last particular result with that one found in the reference ||, we con- 
clude that the intensity-dependent and intensity-independent generalized Jaynes-Cummings 
Hamiltonians have the same eigenstates in the resonant situation. 

b) The Standard Intensity-Dependent Jaynes-Cummings Limit 

The second important limit corresponds to the standard intensity-dependent Jaynes- 
Cummings case, related with the harmonic oscillator system. In this limit we have that 
T = T* — ► 1, 5_ — ► a, B + — > a\ A = uj — uj and £ m +i = (m+ l)hu. Using these 
conditions in Eqs. (CTl) , (^22]) and Eqs. (p|) we can promptly conclude thai; 



#W = (m + l)hu ± hJa 2 LU 2 (m + l) 2 + (u - uo ) 



and 



Im m m 



(2.27) 



(2.28) 



where 



5 



I III 
S m = 



1 + 8^ T $m ■ 



UJ — UJ n 



au(m + 1, 



(2.29a) 
(2.29b) 



Therefore the standard intensity-dependent Jaynes-Cummings eigenstate, written in a ma- 
trix form, is given by 



i 



"1 " 




m) 






. m + 1} . 



m = 0, 1,2,- ■■ . 



(2.30) 



III. TIME EVOLUTION OF THE SYSTEM 



To resolve the the time-dependent Schrodinger equation for intensity-dependent and 
non- resonant Jaynes-Cummings systems 



ih^\*{t)) = (H + H int ) \*{t)) 



we can write the state as 



# (t)> = exp (-m t/n) | %(t)) 



(3.1) 



(3.2) 



and, by substituting this into Schrodinger equation and taking into account the commutation 
property between H D and H int , we obtain 



(3.3) 



Now, we can introduce the evolution matrix Uj(t,0), related with the interaction Hamilto- 
nian, by 



with 



= Ui(t,0) | ^(0)> 



ih^-Vi(t,0) =H int U i (t,0) 
at 



(3.4) 



(3.5) 



that is, in matrix form, written as 

ih 



TV IV 

u 21 u 22 



a 



(3 TB- yjHt 

H 1 B+f t -p 



U21 U22 



(3.6) 



where the primes denote the time derivative. One fast way to diagonalize the evolution 
matrix differential equation is by differentiating Eq. (|3.5| ) with respect to time. After that, 
if we use again the same Eq. fl3.5|), we find 
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Hi„t^Ui(t,0) 
ot 



which can be written as 



u ll ^12 
u 21 u 22 



Cj x 
u 2 



C/„ f/ 12 

^21 ^22 



where 



fajj = a\J(TB_B + nf + (3 2 = \ja 2 # 2 2 + (^A) 



fiwa = a<J{B+BJ)* + (3 2 = J a 2 H 2 + {HA) 2 . 



(3.7) 

(3.8) 

(3.9a) 
(3.9b) 



Now, since by initial conditions Uj(0, 0) = I, then we can write the solution of the evolution 
matrix differential equation ( |3.7|) as 



U,,(t,0) 



cos (&it) sin (uit) C 
sin (o>2^) D cos (^2*) 



(3.10) 



and the C and D operators can be determined by the unitary transformation conditions 

Ut(t, 0) Ufa 0) = U,(t, 0) U](t, 0) = i . (3.11) 

Following the same steps used in the appendix A of the reference we can conclude that 
to satisfy the unitary conditions ( |3 . 1 1| ) these operators must have the form 



C 

b 



H. 



1/4 



B/T^ 



H. 



1/4 



(3.12a) 
(3.12b) 



Therefore, we can write the final expression of the time evolution matrix of the system as 



Ui(t,0) 



cos (uit) sin (obit) C 
— sin (u^i) & cos (o>2^) 



(3.13) 



For Jaynes-Cummings systems an important physical quantity to see how the system 
under consideration evolves in time is the population inversion factor [P,|l|,[l4]] , defined by 



W(t) = a+(t) <7_(t) - <7_(t) &+{t) = Ht) 



(3.14) 



where the time dependence of the operators is related with the Heisenberg picture. In this 
case, the time evolution of the population inversion factor will be given by 



^ = lul(t,0) [* 8) H] U^O) 
at in 1 J 



(3.15) 



and since we have 



a 



&3, Si 



-2a Si <3" 3 . 



(3.16) 
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then Eq. ( [3.15] ) can be written as 



-^- = - r m<t) 



(3.17) 



We can obtain a differential equation with constant coefficients for <J 3 (t) by taking the time 
derivative of Eq. (|3.17|) 



d 2 a 3 {t) 2ia \ dS^t) 



dt 2 



n 



dt 



& 3 {t) + Si(t) 



dt 



(3.18) 



Having in mind that 



and, 



we can conclude that 



Sj, H 



Ui(*,0) 



Si,H = a(3 Si, 6- 3 =2a/3Si<T 3 . 



Si(t) a 3 {t) 



dt h 

Now using Eqs. ( 3-17Q and (|3.21|) into Eq. (|3.18|) we obtain 

d 2 a 3 (t) 



dt 2 



+ & 2 a 3 (t)=F(t) 



where 



e 2 = ^sj 

h 2 

F(t) = ^TJl(t,0) S<U,(f,0). 



(3.19) 
(3.20) 
(3.21) 

(3.22) 

(3.23a) 
(3.23b) 



The Eq. ( p.21|) corresponds to a non- homogeneous linear differential equation for a 3 (t) 
with constant coefficients since Sf and H commute and, therefore, is a constant of the 
motion. The general solution of this differential equation can be written as 

a 3 (t)=a H (t) + a p (t), (3.24) 

and each matrix element of the homogeneous solution, satisfies the differential equation 



+ v 2 a? k (t) = 0, j, fc = l,or2 



dt 2 



with 



hi>i = 2aTB-B + f ] = 2a H 2 
hv 2 = 2aB + B_ = 2aH l . 



(3.25) 



(3.26a) 
(3.26b) 
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The solution of Eq. ( |3.25| ) is given by 

ofiSt) = Vjtt) Cjk + Zj(t) djk , (3.27) 

where 

yj(t) = cos (i>jt) (3.28a) 

Zj(t) = sin(z/,t) , (3.28b) 

and the coefficients Cjk and dj^ can be determined by the initial conditions. 

The matrix elements of the particular solution of the 0*3 (t) differential equation needs to 
satisfy 



(it 1 ' ' • ' 



J jk\ v ) ~2 ±.P 



+ 9 af k {t) = F jk (t) , j, k = 1, or 2 , (3.29) 



and can be obtained by the variation of parameter or by Green function methods, giving 

&m = v- 1 fat) [^m)FAO-m f de^-co , (3.30) 



where we used that the Wronskian of the system of solutions y~j(t) and %(t) is given by i>j. 

After we determine the elements of the F(t)-matrix, it is necessary to resolve the integrals 
in Eq. ( |3.30| ) to obtain the explicit expression of the particular solution. In the appendix we 
show that, using Eqs. (|2.2|), ( [3.13j ), and (|3.23|) , it is possible to conclude that these matrix 
elements can be written as 



ofi(t) = 1 \lTB- {z 2 {t) Qcs(t; z> 2 ,^2,^i) -hif) Gssfa h,<^2,<^i)} H 2 /4 



+ O-v^ET fat) G { s ^{t] *>i,a>iA) - Vi(t) G { cc^ ^lAA)} >fS+fr, (3.31a) 



2 

*f 2 (*) = 1"iWt b - fat) Occ(t; ^2,^2,^1) - m(t) G { sc(ti h^2^i)} \Jh 2 tb 



2 

+ ^rW 4 fat) g ( ss\t; &i,£>i,£>2) + Vi(t) G ( cs\t; OxMM}^ , (3.31b) 



2 

= 1 -y 2 1 \[B~4^W 2 {z 1 {t) g$}(t; KMM) - Vi(t) g { s$(t; u> 2 )} ^B~f^ 

+ Z^Hx* {z 2 (t) Gss(t; z>2, &2, wi) - y 2 (t) G { cs(t; h, ^2, £1)} H 2 /4 , (3.31c) 
= t 1 -o 2 1 \[B~m 2 {z l {t) G^(t; - Vi(t) Gss\t; &i,&i,&2)}h} /4 ' 

+ ih^H 1 ^ fat) g { iJ(t; i>2,ou2,oui) + fa{t) g [ c^(t; z> 2 , c£; 2 , c^)} \Jh 2 TB_, (3.31d) 

where 7 = 4a 2 (3 /h 2 , and the auxiliary functions are given by 

g^(t-p,q,r) = ^ XY (t-p-q,r)±^ XY (t-p + q,r), X, Y = C or S , (3.32) 

with 
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J~cc(t',x,w) = / d£ cos cos(w£) 
Jo 



x 2m w 2n 



i2m+2n+l 

) (2m)!(2n)!(2m + 2n + l) 

ft 



(3.33a) 



J~cs(t] x,w) = / d£, cos(££) sin(tu£) 
Jo 



E (-1: 

m,n=0 
ft 



m+n 



x 2m w 2n+1 



j2m+2n+2 



(2m)! (2n + 1)! (2m + 2n + 2) 

x,w)= f d£ sin (sc£) cos (wg) 
Jo 



(3.33b) 



E (-1: 



m+n 



2m+2n+2 



m,n=0 

J~ss(t', x,w) = / d£ sin sin (w£) 



(2m + 1)! (2n)! (2m + 2n + 2) 



(3.33c) 



E (-1: 



m+n 



^2m,+1^2n+l 



2m+2n+3 



m ,n=o (2m+l)!(2n + l)!(2m + 2n + 3) 

With these results for the particular solution we can conclude that 

^ _ n _ 



(3.33d) 



<tf(0) = 



tit 



Now, using Eqs. (|3.17|) , (|3.24|) , (|3.27|) and the initial conditions, we have 

[*s(0)L 



dt 



2ia r 
T 



Si(0)ff 3 (0) 



(3.34) 

(3.35a) 
(3.35b) 



Therefore, the final expression for the elements of the population inversion matrix of the 
system can be written as 



2ia 



(3.36) 



[£ 3 (t)] y = cos {pit) MO)]^ + — sin (Oit) z/r 1 [S(0)i & 3 (0)\ .. + ££(t) . 

Again, using these final results we can verify two limiting cases, 
a) The Intensity-Dependent Resonant Limit 

The first one corresponds to the intensity-dependent resonant (A = 0). Eqs. (|3.9| ), Q3.13|) , 
(p.26 ) and ( |3.31|) allow us to conclude that, in this case, the evolution matrix of the system 
is given by 



— sin 



Ui(*,0) = 

and the elements of the population inversion of the system is 



cos (l^i tj sin (f £itj C 
Qz> 2 t) & cos ( |z> 2 t 



(3.37) 
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2ia 



Mt)}^ = cos (pit) [a- 3 (0)] y + — sin^t) vr v [S i (0)a 3 (0) 



(3.38) 



b) The Standard Intensity-Dependent Jaynes-Cummings Limit 



This second important limit corresponds to the case of the harmonic oscillator system, 
and in this limit we have that T = — > 1, B- — > a, B + — ► a* and [d,a^] = tkv. 
With these conditions the operators Cj% and Cj 2 commute, and this fact permits to evaluate 
the integrals related with the particular solution of the population inversion elements using 
trigonometric product relations. Using that and the expressions obtained in the appendix, 
after a considerable amount of algebra and trigonometric product relations we can show that 
is possible to write the expressions for the <7^(t)-matrix elements as 



■if- 
.7 



3/4 



i-v x 1 Va~{JC s (t; u 2 ,u x , z> 2 ) - JC s {t; a) 2 , —Ui, z> 2 )} (aa? 

3/4 



%-v x 1 (aa f ) {JC s (t; u 2 , G) x , z>i) - K. s (t; u 2l -Cj x , z>i)} v 7 ^ 7 



aa'a 



-z^ 1 (aa f ) 1 {Ka{t\ u 2 , u>i, i>i) - JC c (t; u 2 , -u h z>i)} (tfa) 



^2i(t) = -v 2 1 \/alaltf {K. c (t;uJ 2 ,u}i,h) +JCc(t]tb 2 ,-ui,u 1 )}v¥ 



7.-1 



a'a 



1/4 



2 

r 

I - 

( 

z— z/ 2 ' " " ' 



{K c (t; lj 2 , tbi, v 2 ) - Ko(t] u 2 , -lu x , z> 2 )} (aa f ) 



1/4 



3/4 



1/4 



= 1 "^^{^5'(*; w 2 , Wi, />i) + /C 5 (t;c^ 2 ,-Wi,z> 1 )} (a f a) 

fc^a") / {/Cs(t; a> 2 , ^2) + fcs(t; & 2 , — u x , v 2 )} V hat a 

where, now, the auxiliary functions are given by 

f sin [(p + q) t] — (p + g) sin (ft) 



£<?(£; p,g,f) 



JC c {t;p,q,f) 



f 2 — (p + (?) 
f cos [(p + g) t] — f cos (ft) 
f 2 — (p + g) 2 



(3.39a) 
(3.39b) 
(3.39c) 
(3.39d) 

(3.40a) 
(3.40b) 



Considering the expressions above we may easily verify that the particular solution for the 
population inversion factor must still satisfy the initial conditions ( |3.34|) . Therefore, in this 
case the final expression of the population inversion factor has the same form given by 
Eq. flQgP , with 



Tii>\ = 2a da 



huji = aJ (aat) 2 + /3 2 



T1O2 = 2a a'a , 

hu 2 = av /(ata) 2 + /5 2 



(3.41a) 
(3.41b) 
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IV. CONCLUSIONS 



In this article we introduced a class of shape-invariant bound-state problems which rep- 
resent two-level systems. The corresponding coupled-channel Hamiltonians generalize the 
intensity-dependent and non-resonant Jaynes-Cummings Hamiltonian. These models are 
not only interesting on their own account. Being exactly solvable coupled-channels prob- 
lems they may help to assess the validity and accuracy of various approximate approaches 
to the coupled-channel problems JTH . 
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Appendix 



In this appendix we show the steps necessary to obtain the explicit expression of the 
elements of the population inversion particular solution. To resolve the integrals in the 
Eq. ( |3.30| ), first we need to determine the elements of the F(t)-matrix. To do that we can 
use Eqs. ( ^23|) , and ( PP to conclude that 



F 12 (t) 



F 2 i{t) 



F 22 {t) 



—7 {cos {(bit) TB_ sin {(b 2 t) & + C sin {(b 2 t) B + T^ cos 

T£>_ cos {(b 2 t) sin {(bit) H^ 4 — H^ 4 sin {(bit) cos {(b 2 t) \J B + T^ 



7 |cos {(bit) TB_ cos {(b 2 t) — C sin (a> 2 ^) B + T^ sin {(bit) Cj 

TB_ cos {(b 2 t) cos {(bit) \JtB_ + H^ 4 sin {(bit) sin {(b 2 t) H^ 4 



7 



7 |cos {(b 2 t) B + T^ cos {(bit) — & sin {(bit) TB^ sin {(b 2 t) 
7 



B+f t cos {(bit) cos {(b 2 t) yB + f^ + H{ 14 sin {(b 2 t) sin {(bit) H 2 1/4 
7 sin {(bit) TB_ cos {(b 2 t) + cos {(b 2 t) B + T^ sin {(bit) C*| 



y\J B + T^ cos {(bit) sin {(b 2 t) H^ 4 — H^ 4 sin {(b 2 t) cos {(bit) \Jtb\ 



(la) 



(lb) 



;ic) 



(Id) 
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where 7 = 4a 2 /3//i 2 , and we used the properties 



C& = C ] C = 1 
C sin {Cj 2 i) = sin (u\t) C 
cos (uit) = cos (u> 2 £) Cfl 



proved in the appendix A of the Ref. [ED], together with the operators relations 



(2a) 
(2b) 
(2c) 

(2d) 

(2e) 



C\B^ 



TB_ C ] = iH 2 1/4 
S+ft C = -C ] \[fB~L = iH\ /A . 



(3a) 
(3b) 

At this point, if we remember that [£j,cDj] = 0, (j = 1, or 2), then we conclude that we 
can use the trigonometric relations involving the product of trigonometric function with argu- 
ments i>jt and Cjjt because, in this case, we know that exp {pjt) exp (±a)jt) = exp [{pj ± Cjj)t]. 
Now, using this fact, the commutators 



z>i,if 2 



z> 2 , Hx = U)2,Hi = , 



(4) 



and the properties (Q), we can show that 



■ 7 



Vi(t) F u (t) = i-\JTB- {cos [(i> 2 - iu 2 )t] sin (ibit) + cos [(z> 2 + u) 2 )t] sin (a>it) } if 2 



1/4 



+ i^-i^ 4 {sin [(z>! — cos (<u> 2 t) — sin [(pi + £>i)t] cos (u) 2 t)} y £? + Tt 



(5a) 



= i^\jTB_ {cos [(z> 2 - & 2 )£] cos (cuit) + cos [(z> 2 + u) 2 )t] cos (uit)} \JfB_ 
+ ^■-^2 / ' 4 {si 11 [(^1 + ^i)^] si n (^2^) — sin \{p\ — Cj\)t\ sin (o) 2 £)} i^ 4 



(5b) 



&(*) F 2 i(t) = ^B+fi {cos [(Pi - £i)t] cos (w 2 t) + cos [{v x + u>i)t] cos {u 2 t)} \fB + fi 



+ ^H^ 4 {sin [(z> 2 + o) 2 )t] sin (tDji) — sin [(z> 2 — a) 2 )i] sin i^ 2 ^ 4 
.7 



(5c) 



y 2 (t) F 22 (t) = i-y B + T^ {cos [(z>i - cDx)*] sin (a) 2 t) + cos [(z>i + &x)t] sin (cj 2 t)} ^ 



1/4 



+ i^-H\^ A {sin [(i> 2 — a) 2 )t] cos (a>it) — sin [(z> 2 + a) 2 )t] cos (u>it)} \JtB^ . (5d) 



In the same way, we can show that 



• 7 



zi(^) F u (t) = i-yTB^ {sin [(P 2 - Cj 2 )t\ sin (u>it) + sin [(z> 2 + u) 2 )t] sin {(bit)} H 2 



1/4 



— z^i^ 4 {cos [(Pi — o>i)t] cos (u) 2 t) — cos [(z>i + cos (a) 2 t)} y 7 £? + Tt (6a) 
-Si (t) A 2 (£) = i^\jTB^ {sin [(z/ 2 - u) 2 )t] cos (u>it) + sin [(z> 2 + u 2 )t] cos (a>it)} \JfB_ 
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— 7^2^ 4 { cos [(At + ^i)t] sin (^2*) — cos — &i)t] sin (a) 2 t)} H\^ a (6b) 
hit) F 21 (t) = -y/B+fi {sin [(z>i - cos (u 2 t) + sin [(^ + uj x )t] cos (w 2 t)} ^fi+ft 

— ^-^i^ 4 {cos [(^2 + &2)t\ sin — cos [(z> 2 — £2)^] sin (u>]t)} (6c) 
z 2 (t) F 22 (t) = i^\J B + f^ {sin [(z>i - u>i)t] sin {u 2 t) + sin [(z>i + tdijt] sin (cu 2 t)} i^ 1/4 

— i-^H^ 4 {cos [(z> 2 — co 2 )t] cos (cDxi) — cos [(z> 2 + iu) 2 )t] cos (£>]£)} \JtB_ . (6d) 

Now, the non-commutativity between the operators C)\ and u 2 imply that to calculate the 
integrals involving the terms given by the Eqs. @ and (H) we need to use the series expansion 
of the trigonometric functions. In this case the integrals can be easily done because the time 
variable can be considered as a parameter factor. Finally, using these results into Eq. Q3.30 ) 
is trivial to find the expression ( |3.31| ) for the matrix elements of the particular solution. 
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